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Mixing-demixing transition in one-dimensional mixtures of fermions and bosons is numerically 
investigated by changing various parameters such as number densities of each component, fermion- 
boson interactions, boson-boson interactions, and hopping energies. In most cases we found clear 
evidences of the mixing-demixing transitions and identified the roles of each microscopic parameter 
in the transitions. Several phase diagrams were obtained. 
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I. INTRODUCTION 

Recent studies of ultracold trapped atomic gases have 
led to the discovery of intriguing physical phenomena and 
have attracted more and more attentions to the atomic 
systems. The trapped atomic gases have some advan- 
tages as a physical system in that any combinations of 
bosons and fermions can be put together on a lattice, mi- 
croscopic parameters are tunable in most cases, and the 
systems can be constructed in various geometries, e.g. 
one-dimensional chain or two-dimensional plane. 

One dimensional systems were constructed with bosons 
0-0 and with bosons and fermions 0. These could be 
the first case that realized one-dimensional boson/boson- 
fermion systems, and triggered theoretical studies 0- 

m 

At a commensurate filling with the lattice periodic- 
ity, the mixtures of spinless fermions and bosons in one- 
dimensional optical lattice were predicted to undergo 
Mott transition at sufficiently large fermion-boson and 
boson-boson interactions 0, • 

In the boson-fermion mixtures at an incommensurate 
filling, on the other hand, the only possible instability 
is mixing-demixing transition. When the repulsive inter- 
actions between fermions and bosons are small, the two 
components are well mixed and the system is in the mix- 
ing state. When the interactions are large enough, each 
component stays away from the other and the system is in 
the demixing state. Cazalilla and Ho studied the mixing- 
demixing transition based on a Luttinger liquid formal- 
ism with a large difference in the number densities of the 
two components |l3j . They derived a criterion for the 
demixing of hardcore bosons and spinless fermions in the 
form, y/vfVFb < |<7/&|/^ 7r where g/b is a fermion-boson 
interaction, Vf is the Fermi velocity, and Vpb = hirpb/Mb 
with Mb being the atomic mass of boson and pb being 
the boson density. Also there is a rigorous proof that 
the system always remains in the mixing state as far as 



the fermion-boson interaction has the same magnitude as 
the boson-boson interaction [lOj . For a slightly different 
but more generalized model, Batchelor et al. found no 
demixing phase in the Bethe Anzatz solution .11.]. 

Under these circumstances we need more evidences for 
the existence of the mixing-demixing transition to draw 
clear phase diagrams. In this paper, we investigated nu- 
merically the mixing-demixing transitions of incommen- 
surate fermion-boson mixtures to complete the phase di- 
agrams, and identified the roles of various microscopic 
parameters in the transition. 

The paper is organized as follows. In Chapt. 2 we 
will introduce the model Hamiltonian and a correlation 
function we used as an indicator to determine the phases. 
Chapter 3 is devoted to the description of calculation 
results and Chapt. 4 to the conclusion of the present 
study. 



II. MODEL 

We consider a mixture of Nf spinless fermions and Nb 
bosons on an iV-site lattice. We employed Bose-Fermi 
Hubbard Hamiltonian to describe the interacting parti- 
cles, which is defined by 



U = ^2 \t a {a i a , i a a ,i+i + h.c.)+ p a n a 

i a—f.b 

~u bb 



E 



-nb,i(n b ,i - 1) + UfbUf^nb^ 



, (1) 



where a a i and a a .i are respectively creation and annihi- 
lation operators for fermions (a — /) or bosons (a = b) 
on the i-th site, and n a ^ — j-a a) j. Hopping energy and 
chemical potential are denoted respectively by t a and p, a . 
Ubb and Ufb are on-site boson-boson and fermion-boson 
repulsive interactions respectively. 
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To observe the mixing-demixing transition, we define 
a correlation function 

C = ({n b ,i-i+n bA+1 )n b . i ) - ((n f . i _ 1 +n f . l + n fA+1 )n b . i ). 

(2) 

The first term gives the correlation in neighboring 
bosons, and the second term expresses the correlation 
between bosons and fcrmions on the same and neighbor- 
ing sites. In the mixing states fermions and bosons are 
spatially mixed and move throughout the system, while 
in the demixing states fermions and bosons reside sepa- 
rately and move in each restricted area. The correlation 
function C changes drastically at the transition, reflect- 
ing these characteristics. It tends to be smaller, or neg- 
atively larger, in mixing states than in demixing states. 
The reason we omitted (n b ^n b ^) from the first term is 
that this sometimes dominates the correlation function 
when multiple bosons stay on the same site, mostly in the 
demixing states, and hides subtle behavior of the other 
terms. Many physical quantities change at the transition 
and we actually calculated other various physical quanti- 
ties including the i — j dependence of (n^n/j). However 
we found the correlation function C exhibited the most 
drastic change near the transition point and decided to 
use it as an indicator of the transition. 




FIG. 1: The snapshot of fermions and bosons with Sp — 0, 
tf = 1 and Ubb = 0. In (a) we set Ufb = 0.5 and in (b) 
Ufb ~ 6. The upper graphs in (a) and (b) show the snapshots 
of bosons and the lower graphs show those of fermions. The 
real space is indicated by X axis, r denotes the imaginary 
time axis, and nj, (n/) axis presents the number of bosons 
(fermions) on each site. The variation of the darkness is just 
a guide for the eyes. 



III. CALCULATION AND RESULT 

We performed Monte Carlo simulations with the world- 
line algorithm [T3.[T5| on a one-dimensional iV-site lattice 
with periodic boundary conditions. We changed in the 
simulations the number of fermions and bosons, Nf and 
N b respectively, while fixing the total number of particles 
Nf + N b to 14. We used the hopping energy t b of bosons 
as an energy unit. Temperature was fixed to T = 0.08 
and the Trotter decomposition number L to 100. The 
number of sites N was set to 30. 

In the following two subsections we will show the 
mixing-demixing transitions by changing various param- 
eters with several choices of the difference in the number 
densities of the fermions and bosons, Sp = (Nf — N b )/N. 

A. Changing interactions 

At first we observed the mixing-demixing transitions 
by changing the interactions with the hopping energy ra- 
tio fixed to tf/tb = 1. Figure 1 presents typical snap- 
shots of the fermions and bosons on the (real space)- 
(imaginary time) plane with U bb = and Sp = 0. As far 
as the interaction Uf b between the fermions and bosons 
is small, both components are well mixed as in Fig. 1(a), 
while they are separated at a large Uf b as in Fig. 1(b). 
Due to the absence of the boson-boson interaction, more 
than one boson can occupy the same site with no ad- 
ditional energy cost and hence the bosons can provide 
a sufficiently large space for the fermions. Next we ob- 



served the mixing-demixing transitions in more quanti- 
tative way. Figure 2 shows the correlation function C 
as a function of the fermion-boson interaction Uf b with 
various dp and U bb . U bb was set to in Fig. 2(a), 0.4 
in (b) and 1 in (c). The difference in the number densi- 
ties of the fermions and bosons, Sp = (Nt — N b )/N, was 
chosen to be 0, 2/30, 4/30, 6/30, and 8/30. As we can 
see in the figure, there is a characteristic uprise in ev- 
ery curve, indicating the transitions from the mixing to 
the demixing states. We repeated the same calculations 
in different system sizes and found the uprising behav- 
ior became sharper consistently in the larger size. The 
transition curves were still broad to locate the transition 
point accurately, although we can obtain from the fig- 
ure an evidence for the existence of the mixing-demixing 
transitions and determine the approximate location of 
the mixing and demixing phases in a parameter space. 
Finite size scaling should be done in future to identify 
the precise transition point. 

Several features of Fig. 2 should be noted. Firstly the 
transition occurs at any Sp, not just at relatively large 
Sp as indicated in [l3j . It occurs even when Sp = 0. 
Secondly the transition point seems to shift to larger U f b 
side as the interboson interaction U bb increases. 

The second feature can be understood in the following 
way. In the demixing state, the bosons make a relatively 
small cluster by multiple occupation of sites as far as U bb 
is small, leaving a sufficiently large space for fermions to 
move around. Suppose we increase U bb in this state. It 
becomes hard for the bosons to stay in the small cluster. 
The boson cluster hence spreads to a certain extent to 
reduce the energy cost due to U bb , which restricts the 
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FIG. 2: The correlation function C as a function of Ufb with 
Ubb — in (a), 0.4 in (b) and 1 in (c). tf is set to 1. 
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FIG. 3: Current-current correlation function J of fermions 
(a) and bosons (b) as a function of Ufb with Ubb = and 
tf = 1. The function J is normalized by J (Ufb = 0). 
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FIG. 4: Phase diagram in Nf — Nb versus Ufb with Ubb = 
(a) and U bb = 1 (b). tf is set to 1. "A" and "□" present the 
mixing and the demixing states, respectively. "*" is marked 
for undetermined phase points. 



where j a presents the current density of the fermions 
when a = f and that of the bosons when a = b. The 
correlation function J of each kind of quantum particles 
is expected to be large in the mixing states where bosons 
and fermions extend over the system, and to be small 
in the demixing states where bosons and fermions are 
confined in their own restricted areas. Fig. |3| shows the 
J as a function of Ufb with tf — 1 and Ubb — 0. Every 
curve in the figure suggests the transition from mixing to 
demixing and the transition point is located somewhere 
in between Ufb = 1 and 3, which is consistent to what 
we have seen in Fig. 2(a). 

Let us now draw the phase diagrams. Figure 4 shows 
the phase diagram in the Sp-Ufb plane. "A" and "□" 
denote the mixing and demixing states respectively that 
can be determined clearly from the data in Figs. 1-3. As 
stated in the above, the transition curves are too broad 
to identify the precise location of the transition point 
and therefore we used the symbol "*" in the phase dia- 
gram to present a gray zone where we cannot determine 
whether the particles are mixed or demixed. However it 
is still clear that increasing Ubb drives the transition area, 
presented by *, toward the larger Ufb side. 



fermions to gain the kinetic energy. In other words, Ubb 
helps mixing, while Ufb helps demixing. Although we 
did not introduced interfermion interactions Uff in the 
model, it is obvious Uff would have the same effect as 
Ubb- 

We also observed the transitions by measuring an- 
other physical quantity, current-current correlation func- 
tion J a =f.b in the zero-frequency limit [l5l llrJ]: 

J a = lim (J a (uj)j a (-uj)), (3) 



B. Changing tf 

We then observed the transitions by changing tf. Fig- 
ure[5]shows the snapshot of the bosons and fermions as in 
Fig. 1, although the hopping energy tf is increased here. 
By comparing (a) and (b), we see that the fermions are 
invading the boson area and making the system mixed 
as tf becomes larger. 

Figure shows the correlation function C versus U fb 
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FIG. 5: The snapshot of fermions and bosons with 5p = 6/30, 
Ufb = 6 and Ubb ~ 0. In (a) we set tf = 2 and in (b) tf = 7. 
See the caption of Fig. 1 for other details of the figure. 

for two different values of tf. Figure El (a) is equivalent 
to Fig. 2(a) with the same parameters, and Fig. El (b) 
shows the correlation function for larger tf. As we see 
in the figure the transition points move to the larger U fb 
side as tf increases. In general, when hopping energy 
increases, particles need more space to gain the kinetic 
energy. This tendency is more drastic in fermions than 
in bosons as far as the on-site interboson interaction is 
finite. With the increase of tf in the demixing phase, the 
fermions try to invade the boson territory. Namely, tf 
helps mixing of the two components. For this reason the 
mixing-demixing transition occurs at a larger U fb when 
tf becomes large as in Fig. 

The same behavior can be found in Fig. [7\ where the 
correlation function J of the fermions (a) and the bosons 
(b) is calculated. Figure should be compared with Fig. 
3. As tf increases the transition points shift to the larger 
Ufb side. 

We also drew the phase diagram in Fig. El with a 
larger tf than used in Fig. 4 (a). In comparison of Fig. El 
with Fig. 4(a) we see the transition area, represented by 
*, clearly shifts to the larger Ufb side. 

IV. CONCLUSION 

We studied the mixing-demixing transitions of 
fermion-boson mixtures in one dimension by changing 
various parameters. Monte Carlo calculation showed the 
mixing-demixing transition occurred at any dp, the dif- 
ference in the number densities of the two components, 
as the repulsive interactions between the fermions and 
bosons increase. The schematic phase diagram is given 
in Fig. E| Also we confirmed the role of each parameter. 
The fermion-boson interaction Ufb helps demixing, while 
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FIG. 6: The correlation function C as a function of Ufb with 
Ubb = 0. tf is 1 in (a) and 4 in (b). (a) is equivalent to Fig. 
2(a). 
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(a) (b) 

FIG. 7: Correlation function J of the fermions (a) and the 
bosons (b) as a function of Ufb with Ubb = and tf = 4. The 
function J is normalized by J (Ufb ~ 0). 

the boson-boson interaction Ubb and the hopping energy 
tf help mixing. The phase boundary of Fig. El therefore 
shifts to the larger Ufb side as Ubb on tf increases. 

In our simulations we could not obtain reliable results 
at large interactions because of the restriction on the 
computational time. We performed the simulations for 
fermions and hardcore bosons, i.e. with Ubb — oo, and 
found no demixing phase as far asUfb < 20. Although we 
have no numerical evidence on the absence of the demix- 
ing phase at Ufb > 20, we believe the fermions and the 
hardcore bosons would not be demixed at any U fb for 
the following reason. Suppose we have demixed fermions 
and hardcore bosons with a finite Ufb- In the vicinity 
of the boundary, the fermions and bosons move around, 
interacting with each other. The fermions (bosons) can 
only move to empty sites or sites occupied by the bosons 
(fermions), which means that it is energetically favored 
for each component of the particles to invade the area 
of the other component rather than to stay in its own 
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area. In other words the demixing state of the fermions 
and the hardcore bosons would be unstable against mix- 
ing. In fact there is a rigorous proof that a mixture of 
fermions and bosons remains in the mixing state when 
Ufb = Ubb 10], indicating finite Ufb cannot demix the 
mixture of fermions and hardcore bosons. 



FIG. 8: Phase diagram in Nf — JVj, versus Ufb with Ubb = 
and tf — 4. "A" and "□" present the mixing and the 
demixing states, respectively. " *" is marked for undetermined 
phase points. 
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FIG. 9: Schematic phase diagram of the mixing and the 
demixing phases. 
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